library (foreign)
library (survival)
library (arm)
library (pscl)
library(apsrtable)


######################################################
#
#   Negative Binomial for MID Initiation
#   Thursday, July 26 1:55 p.m.
#######################################################



TermLimitsData<-  read.dta("/Users/Jeff/Dropbox/TermLimits/D6/Analysis/Data/TermLimitsAnalysisData.dta")

summary(TermLimitsData)
#TermLimitsData <- subset(TermLimitsData, year >= 1960)
#attach(TermLimitsData)



##  Unconditional   ##

mod1<- glm.nb (countinit ~ lameduck + hawkish +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


mod2<- glm.nb (countinit ~ lameduck + hawkish2 +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


mod3<- glm.nb (countinit ~ lameduck + milservice +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


apsrtable(mod1,mod2,mod3)




##  Conditional   ##

mod1<- glm.nb (countinit ~ lameduck + hawkish + lameduck_hawkish +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


mod2<- glm.nb (countinit ~ lameduck + hawkish2 + lameduck_hawkish2 +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


mod3<- glm.nb (countinit ~ lameduck + milservice + lameduck_mil +
                           rivalry + numbord + parliament + gender +
                                   peace + peace2 + peace3,data=TermLimitsData)


apsrtable(mod1,mod2,mod3)



